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Abstract 

In recent years single molecule force spectroscopy has opened a new avenue to provide profiles of the 
complex energy landscape of biomolecules. In this field, quantitative analyses of the data employing 
sound theoretical models, have played a major role in interpreting data and anticipating outcomes 
of experiments. Here, we explain how by using temperature as a variable in mechanical unfolding of 
biomolecules in force spectroscopy, the roughness of the energy landscape can be measured without 
making any assumptions about the underlying reaction coordinate. Estimates of other aspects of the 
energy landscape such as free energy barriers or the transition state (TS) locations could depend on 
the precise model used to analyze the experimental data. We illustrate the inherent difficulties in 
obtaining the transition state location from loading rate or force-dependent unfolding rates. Because 
the transition state moves as the force or the loading rate is varied, which is reminiscent of the Hammond 
effect, it is in general difficult to invert the experimental data. The independence of the TS location 
on force holds good only for brittle or hard biomolecules whereas the TS location changes considerably 
if the molecule is soft or plastic. Finally, we discuss the goodness of the end-to-end distance (or 
pulling) coordinate of the molecule as a surrogate reaction coordinate in a situation like force-induced 
ligand-unbinding from a gated molecular receptor as well as force-quench refolding of an RNA hairpin. 
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I. INTRODUCTION 


The energy landscape, which projects the multidimensional conformational space of biopoly¬ 
mers onto low dimension, has played a critical role in visualizing their folding routes [lH3j. To 
account for the rapid, reversible folding and unfolding of proteins, it is suspected that the fold¬ 
ing energy landscape of many evolved proteins is relatively smooth, which allows for an efficient 
navigation of the landscape. To be more precise, the gradient of the energy landscape A F(x) to¬ 
wards the native basin of attraction (NBA), corresponding to the driving force, is “large” enough 
that during the folding process the biomolecule does not get kinetically trapped in local minima 
(competing basins of attraction (CBA)) for arbitrarily long times. Here, F(y) is expressed as 
a function of a non-unique variable, namely, the structure overlap function y, an order param¬ 
eter that measures how similar a given conformation is to the native state. However, perfectly 
smooth energy landscapes are difficult to realize because of energetic and topological frustration 
m ■ In proteins, the hydrophobic residues prefer to be sequestered in the interior while polar 
and charged residues are better accommodated on the surfaces where they can interact with 
water. Often these conflicting requirements cannot be simultaneously satisfied and hence pro¬ 
teins can be energetically “frustrated”. It is clear from this description that only evolved or well 
designed sequences can minimize energetic frustration. Even if a particular foldable sequence 
minimizes energetic conflicts, it is nearly impossible to eliminate topological frustration which 
arises due to chain connectivity [SEj. Topological frustration refers to the conflict between 
local and global packing of structures. Both sources of frustration, energetic and topological, 
render the energy landscape rugged on length scales that are larger than those associated with 
secondary structures (~ (1 — 2) nm). 

Thus, the free energy, projected along a ID coordinate, is rough on certain length scale 
and may be globally smooth on a larger scale. Under the assumption that the characteristic 
roughness e has a Gaussian distribution (P(e) ~ e -e 2 / 2 e 2 ^ the overall transition time from the 
unfolded basin to NBA over free energy barrier A F* and t(/3) = r 0 e^ Ai ' t may be written as 

poo POO 

tu^ f = t{0) / dee^Pie) « r(/3) / dee^e^^ (1) 

J kgT J kgT 

where r(/3) is the temperature dependent Arrehius-like transition time from unfolded (U) to 
folded state (F) over a single barrier in the perfectly smooth one-dimensional coordinate, and e is 
the average value of ruggedness. The last part of the equation becomes valid at low temperatures. 
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The additional factor e^ 2<;2 / 2 in Eq§ which slows down the folding time, was derived in an 
elegant paper by Zwanzig [8] and was also obtained in [9j ;10] by analyzing the dynamics of 
Derrida’s random energy model m If folding takes place in a rough energy landscape then the 
dependence of the characteristic time scale on length l may be estimated as r = t$s ~ (10 — 100) 
ns when l & (1 — 2) nm where Tss ~ l 2 /D is the diffusion time with a diffusion constant on the 
order of (1CT 7 —10 -6 ) cm 2 /sec, and r = when l « L where L is the effective contour length 
of the biomolecule. Given the crude physical picture, the estimate of tss is not inconsistent with 
the time needed to form a-hclices or /3-hairpin. 

With the possibility of mechanically manipulating biological molecules, one molecule at a 
time, it is becoming possible to probe the features of their energy landscape (such as roughness 
and the transition state location) that are not easy to measure using conventional experiments. 
Such experiments, performed using Laser Optical Tweezers (LOTs) (12[ '13] or Atomic Force 
Microscopy (AFM) [13], have made it possible to mechanically unfold proteins [Tol [2T| . RNA 
[22H2E], and their complexes [22H33], or initiate refolding of proteins [3B] and RNA [251127] . These 
remarkable experiments show how the initial conditions affect refolding and also enable us to 
examine the response of biological molecules over a range of forces and loading rates. In addition, 
fundamental aspects of statistical mechanics, including non-equilibrium work theorems [36, 37], 
can be rigorously tested using single molecule experiments [38., 39j. Here, we are concerned with 
using the data and theoretical models to extract key characteristics of the energy landscape of 
biological systems. 

The crude physical picture of folding in a rough energy landscape (Fig |T|) is not meaningful 
unless the ideas can be validated experimentally, which requires direct measurement of the 
roughness energy scale e, absolute value of the barrier height, etc. In conventional experiments, 
in which folding is triggered by temperature, it is difficult to measure To and A F$ even when 
(3e = 0 [IQ]. We proposed, using theoretical methods, that (3e can be directly measured using 
forced-unfolding of biomolecules and biomolecular complexes. The Hyeon-Thirumalai (HT) 
theory j3Tj showed that if unbinding or unfolding lifetime (or rates) are known as a function of 
the stretching force (/) and temperature (T) then e can be inferred without explicit knowledge of 
Tq or A FK Recently, the loading-rate dependent unbinding times of a protein-protein complex 
using atomic force microscopy (AFM) at various temperatures have been used to obtain an 
estimate of e [32]. Similarly, variations in the forced-unfolding rates as a function of temperature 
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of Dictyostelium discoideum filamin (ddFLN4) were used to estimate e [03]. The variation in 
unbinding or unfolding rates of proteins as a function of / and T provides an opportunity to 
obtain quantitative estimates of the energy landscape characteristics. 

Single molecule force spectroscopy can also be used to measure force-dependent unfolding 
rates from which the location of the transition state (TS) in terms of the spatial extension (R) 
can be computed. This procedure is a not straightforward because, as shown in a number of 
studies [HH3B], the location of the transition state changes as / changes unless the curvature 
of the free energy profile at the TS location is large, i.e, the barrier is sharp. The extent to 
which the TS changes depends on the load. By carefully considering the variations of force 
distributions it is possible to obtain reliable estimates of the TS location |43j . Here, we review 
recent developments in single molecule force spectroscopy that have attempted to obtain the 
energy landscape characteristics of biological molecules m 0EHM]. Using theoretical models 
that consider dynamics in higher dimensions we also point out some of the ambiguities in 
interpreting the experimental data from dynamic force spectroscopy. 


II. THEORETICAL BACKGROUND 

Single molecule mechanical unfolding experiments differ from conventional unfolding exper¬ 
iments in which unfolding (or folding) is triggered by varying temperature or concentration of 
denaturants or ions. In single molecule experiments folding or unfolding can be initiated by pre¬ 
cisely manipulating the initial conditions. In both forced-unfolding and force-quench refolding, 
the initial conformation, characterized by the extension of biomolecule, is precisely known. By 
contrast, the nature of the unfolded states, from which refolding is initiated, is hard to describe 
in ensemble experiments. For RNA and proteins, whose energy landscape is complex |3, 55], 
details of the folding pathways can be directly monitored by probing the time dependent changes 
in the end-to-end distance R(t) of individual molecules. Analysis of such mechanical folding and 
unfolding trajectories allows one to explore regions of the energy landscape that are difficult to 
probe using ensemble experiments. 

Force experiments can measure the extension of the molecule as a function of time. There 
are three modes in which stretching experiments are performed. Most of the initial experiments 
were performed by unfolding biomolecules (especially proteins) by pulling on one of the molecule 
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at a constant velocity while keeping the other end fixed m E21E31 [56]. More recently, it has 
become possible to apply constant force on the molecule of interest using feed-back mechanism 
26 . 3oI [o7l Ia9l . In addition, force-quench experiments have been reported in which the forces are 
decreased or increased linearly [26] 2? : , [35]. It is hoped that a combination of such experiments 
can provide a detailed picture of the complex energy landscape of proteins and RNA. In all the 
modes, the variable conjugate to / is the natural coordinate that describes the progress of the 
reaction of interest (folding, unbinding or catalysis). If there is a energy barrier confining the 
molecular motion to a local minimum, whose height is greater than fcgT, then a sudden increase 
(decrease) of extension (force) signifies the transition of the molecule over the barrier. A rip in 
the force-extension curve (FEC) is the signature of such a transition. Surprisingly, for proteins 
and RNA it has been found that the portions of the FEC between rips can be quantitatively 
fit using the semi-flexible or worm-like chain model pj6], 22] [23], 56] ED] {6J.J . From such fits, the 
global polymeric properties of the biomolecule, such as the contour length and the persistence 
length can be extracted mm- 

Single molecule pulling experiments provide distributions of the unfolding times (or unfold¬ 
ing force) by varying external conditions. The objective is to construct the underlying energy 
landscape from such measurements and from mechanical folding or unfolding trajectories. How¬ 
ever, it is difficult to construct all the features of the energy landscape of biomolecules from 
FEC or mechanical folding trajectories that report only changes at two points. For example, 
although the signature of roughness in the energy landscape may be reflected as fluctuations in 
the dynamical trajectory it is difficult to estimate its value unless multiple pulling experiments 
are performed. We had proposed that the power of single molecules can be more fully realized 
if temperature (T) is also used as an additional variable [TT] 0Sj, 62] [53]. By using T and / 
it is possible to obtain a phase diagram as a function of T and / that can be used to probe 
the nature of collapsed molten globules which are invariably populated but are hard to detect 
in conventional experiments. We also showed theoretically that the roughness energy scale (e) 
can be measured [H] if both / and temperature (T) are varied in single molecule experiments. 
The effect of e manifests itself in the 1/T 2 dependence of the rates of force-induced unbinding 
or unfolding kinetics. In the following subsections we first review the theoretical framework to 
describe the force-induced unfolding kinetics and show how the 1/T 2 dependence emerges when 
the roughness is treated as a perturbation in the underlying energy profile. 
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Bell model: Historically, phenomenological descriptions of the force-induced fracture of ma¬ 
terials and unbinding of adhesive contacts were made by Zhurkov [63] and Bell [42] , respectively, 
long before single molecule experiments were performed. In the context of ligand unbinding from 
a binding pocket, Bell [42] conjectured that the kinetics of bond rupture can be described using 
a modified Eyring rate theory [OS], 


k B T (E t_^ f )/ kBT 
h 


( 2 ) 


where k E is the Boltzmann constant, T is the temperature, h is the Planck constant, and k is the 
transmission coefficient. In the Bell description the activation barrier is reduced by a factor 
7 / when the bond or the biomolecular complex is subject to external force /. The parameter 
7 is a characteristic length of the system under load and specihes the distance at which the 
molecule unfolds or the ligand unbinds. The prefactor is the vibrational frequency of a 
single bond. The Bell model shows that the unbinding rates increase when tension is applied to 
the molecule. Although Bell’s key conjecture, i.e., the reduction of th activation barrier due to 
external force, is physically justified, the assumption that 7 does not depend on the load is in 
general not valid. In addition, because of the multidimensional nature of the energy landscape 
of biomolecules, there are multiple unfolding pathways which require modification of the Bell 
description of forced-unbinding. It is an oversimplification to restrict the molecular response to 
the force merely to a reduction in the free energy barrier. Nevertheless, in the experimentally 
accessible range of loads the Bell model in conjunction with Kramers’ theory of escape from a 
potential well have been remarkably successful in fitting much of the data on forced-unfolding 
of biological molecules. 

Mean first passage times: In order to go beyond the popular Bell model many attempts have 
been made to describe the unbinding process as an escape from a free energy surface in the pres¬ 
ence of force jniEniE 2 iEniEZ!. This is traditionally achieved by a formal procedure that adapts 
the Liouville equation describing the time evolution of the probability density representing the 
molecular configuration in phase space. 

For the problem at hand, one can project the entire dynamics onto a single reaction coordi¬ 
nate provided the relaxation times of other degrees of freedom are shorter than the time scale 
associated with the presumed order parameter of interest [ 68 l [69]. In applications to force- 
spectroscopy, we assume that the variable conjugate to / is a reasonable approximation to the 
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reaction coordinate. The probability density of the molecular configuration, p(x,t |xo), whose 
configuration is represented by order parameter x at time t, obeys the Fokker-Planck equation. 


a - £ ^r A = c FP (*M x ,tM = (J) + 1 PM* o). 


( 3 ) 


where D(x) is the position-dependent diffusion coefficient, and F(x) is an effective one¬ 
dimensional free energy, x 0 is the position at time t — 0. If the initial distribution is given 
by p(x 0 ,t = 0|x 0 ) = S(x — x 0 ) the formal solution of the above equation reads p(x,t\x 0 ) = 
e tCFP 8(x — x 0 ). If we use an absorbing boundary condition at a suitably defined location, the 
probability that the molecule remains bound (survival probability) at time t is 


S(xo,t) = I dxp(x, t\xo) = I dxe tCpp S(x — xq). 


( 4 ) 


In terms of the first passage time distribution, Pfp(xo, t) = —dS(x o, t)/dt, the mean first passage 
time can be computed using, 


t(x 0 ) 


dt [tp FP (x 0 ,t)] 


poo 

/ dt 

dS(x 0 ,t ) 

Jo 

dt 


dt J dxe tCFP ^5(x — xq) 
dt [ dxS(x — x^e^r^ 


( 5 ) 


where C^ FP is the adjoint operator. In obtaining the above equation we used S(xo,t = oo) = 0 
and integrated by parts in going from the second to third line. By operating £p P (.x' 0 ) on both 
sides of Eq|5] and exchanging the variable x with x 0 we obtain 


C j FP (x)r(x) = e F ^ x)/kBT ^-D(x)e- F ^ )/kBT ^-r(x) = -1. 


( 6 ) 


The rate process with reflecting boundary d x r{a ) = 0 and absorbing boundary condition r(b) = 0 
in the interval a < x < b, leads to the well-known expression for mean first passage time, 

t(x) = [ dye F{v)/kBT [ V dze~ F{z),kBT . (7) 

Jx D (y) Ja 

Diffusion in a rough potential: In the above analysis the ID free energy profile F(x) that 
approximately describes the unfolding or unbinding event is arbitrary. In order to explicitly 
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examine the role of the energy landscape ruggedness we follow Zwanzig and decompose F(x) 
into F(x) = F 0 (x)+Fi(x) [8]. where Fq(x) is a smooth potential that determines the global shape 
of the energy landscape, and fq(x) is periodic ruggedness superimposed on F 0 (x). By taking the 
spatial average over F^x) using (e ±/3Fl ^); = j dxe^ Fl ^ x \ where l is the ruggedness length 
scale, the associated mean first passage time can be written in terms of the effective diffusion 
constant D*(x ) as, 


D*(x) = 


t(x) 


D(x) 


( e /m(x)^ e -/3Fi(x)y 
[ dye Fo{y)/kBT 


D*{y) 


dze ~ F ^ z)/kBT . 


( 8 ) 


An inversion of roughness potential, i.e., Fi —F\ does not alter D*(x). In the presence of 
roughness D*(x ) < D(x). When /3 2 (Ff(x)) = (3 2 e 2 is small the effective diffusion coefficient 
can be approximated as Dq ps .D 0 exp (—/3 2 e 2 ) where D 0 is the bare diffusion constant. If P{F \) 
is a Gaussian then this expression is exact. The coefficient associated with 1/T 2 behavior is 
due to the energy landscape roughness provided the extension is a good reaction coordinate. 
The dependence of e~ t2 ‘! k B T ~ in D* suggests that the diffusion in a rough potential can be 
substantially slowed even when the scale of roughness is not too large. 

Barrier crossing dynamics in a tilted potential: In writing the one-dimensional Fokker-Planck 
equation (Eqj3]) we assumed that the order parameter £ is a slowly changing variable. This 
assumption is valid if the molecular extension, in the presence of /, describes accurately the 
conformational changes in the biomolecule. 

Following the Bell’s conjecture we can replace F(x) by F(x) — fx in which / ’’tilts” the free 
energy surface. Thus, in the presence of mechanical force Eq|8] becomes 


k (/) — r (/; x ) 


dye :( F °( y ) 


-fy)/k B T_ 


D*(y) 


dze~ {Fo{z) ~ fz)/kBT . 


(9) 


As long as the energy barrier is large enough (see Fig|l]) Eq|9] can be further simplified using 
the saddle point approximation. The Taylor expansions of the free energy potential Fo(a;) — f x 
at the barrier top and the minimum result in Kramers’ equation mm, 


k (/) = r(f) 


27 tIcbT 


0 /9(AF|(/)-/Axt(/)) 


D*mu b (f)u t s{f ) 

( e «AF*(/)—,A*>(,)) (10) 

\UJb{J >UJts{J ) J 







where ojb and c o ts are the curvatures of the potential, |^F 0 (a;)|, at x = Xb and x ts , respectively, 
the free energy barrier AFq = F 0 (x ts ) — F 0 (xb), m is the effective mass of the biomolecule, / is 
the friction coefficient, and Ax* = x ts — %b- 

In the presence of /, the positions of the transition state x ts and bound state Xb change because 
unbinding kinetics should be determined using F 0 (x) — fx and not F 0 (x) alone. Because x ts 
and Xb satisfy the force dependent condition Fq(x) — / = 0, it follows that all the parameters, 
A£*(/), uu(f), and u b (f), are intrinsically /-dependent. Depending on the shape of the free 
energy potential F 0 (x), the degree of force-dependence of Ax*, u> ts , and u>b can vary greatly. 
Previous theoretical studies |H EH E21 E2] have examined some of the consequences of the 
moving transition state. In addition, simulational studies [441 146] in which the free energy 
profiles were explicitly computed from thermodynamic considerations alone clearly showed the 
change of Ax * when / is varied. These authors also provided a structural basis for transition 
state movements in the case of unbinding of simple RNA hairpins. The nontrivial coupling of 
force and free energy profile makes it difficult to unambiguously extract free energy profiles from 
experimental data. In order to circumvent some of the problems, Schlierf and Ricf have used 
Eq|9] to analyze the load-dependent experimental data on unfolding of ddFLN4 and extracted 
an effective one dimensional free energy surface F(x) without making additional assumptions. 
The results showed that the effective free energy profile is highly anharmonic near the transition 
state region m 

Forced-unfolding dynamics at constant loading rate: Many single molecule experiments 
are conducted by ramping the force over time [M, [16] [23] 08]. In this mode the load on the 
molecule or the complex increases with time. When the force increases beyond a threshold value, 
unbinding or bond-rupture occurs. Because of thermal fluctuations the unbinding events are 
stochastic and as a consequence one has to contend with the distribution of unbinding forces. 
The time-dependent nature of the force makes the barrier crossing rate also dependent on t. For 
a single barrier crossing event with a time-dependent rate k(t), the probability of the barrier 
crossing event being observed at time t is P[t) = k(t)S(t ), where the survival probability S(t) 
that the molecule remains folded is given as S(t) = exp J* drk(r )j. 

When the molecule or complex is pulled at a constant loading rate (r/) the distribution 
(P(f)) of unfolding forces is asymmetric. The most probable r/-dependent unfolding force (/*) 
is often used to determine the TS location of the underlying energy landscape with the tacit 
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assumption that the TS is stationary. When rf — df/dt is constant, the probability of observing 
an unfolding event at force / is written as, 


P(f) — k(f) exp 
r f 


f f 1 
/ df'-k(f) 

/o r f 


( 11 ) 


The most probable unfolding force is obtained from dP(f)/df\f=f* = 0, which leads to 


r = ksT nog ( _ r f Axt (n _ 

A xHf*) \v D (f*)e~d AF ^k B T 


( 


log 


1 + /’ 


A x*(f* 

+ log(e' SF ') i (e-' 5F '),}, 


(Axl)'(f) _ (W) (/*) + vtfj*) 


\ 




A x*(f* 


Mf*) A x*(f* 


( 12 ) 


where AFq = F 0 (x ts (f)) — Fo(xo{f)), 1 denotes differentiation with respect to the argument, 
A x(f) = Xts(f) — xq (/) is the distance between the transition state and the native state, and 
z/d(/) = u 0 (f)uts(f)/ 2717 . Note that AFq, and Ax$ depend on the value of / [3T, 05J H 6 j . 
Because f* changes with r/, Ax 1 obtained from the data analysis should correspond to a value at 
a certain /*, not a value that is extrapolated to f* = 0. Indeed, the pronounced curvature in the 
plot of f* as a function of logrj makes it difficult to obtain the characteristics of the underlying 
energy landscape using data from dynamic force-spectroscopy without a reliable theory or a 
model. If AFg, and Ax* are relatively insensitive to variations in force, the second term 


on the right-hand side of Eq 12 would vanish, leading to f* oc (k B T/Ax*) logrq- [50j. If the 
loading rate, however, spans a wide range so that the force-dependence of A Fq, v d , and Aa;* 
are manifested, then the resulting f* can substantially deviate from the linear dependence on 
log rf. Indeed, it has been shown that for a molecule or a complex known to have a single free 
energy barrier, the most probable rupture force f* obeys (1 — f*/f c ) ~ (log r f) v with f c being a 
critical force in the absence of force. The effective exponent v should lie in the range 0.5 < v < 1 
m ■ And a fit with v < 0.5 implies that the unbinding dynamics cannot be explained with a 
one-dimensional barrier picture of crossing d. The precise value of v depends on the nature 
of the underlying potential and is best treated as an adjustable parameter. 
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III. MEASUREMENT OF ENERGY LANDSCAPE ROUGHNESS 


In the presence of roughness we expect that the unfolding kinetics deviates substantially 
from an Arrhenius behavior. By either assuming a Gaussian distribution of the roughness 
contribution, P(Fi) oc e ~ F i / 2e2 ; or simply assuming (3F\ <C 1 and (Fi) = 0, (Ff) = e 2 , one can 
further simplify EqjlOjto 

log*(/)/*„ = -(A F* - f ■ A x*)/k B T - e 2 /k 2 B T 2 . (13) 


This relationship suggests that the roughness scale e can be extracted if logfc(/) is measured 
over a range of temperatures. Variations in temperature also result in changes in the viscosity, r), 
and because 1 oc rj, corrections arising from the temperature-dependence of r) have to be taken 
into account in interpreting the experiments. It is known that rj for water varies as exp (A/T) 
over the experimentally relevant temperature range (5°U < T < 50°U) [75]. Thus, we expect 
log k(f , T) = a + b/T — e 2 /T 2 . The coefficient of the 1/T 2 term can be quantified by performing 
force-clamp experiments at several values of constant temperatures. In addition, the robustness 
of the HT theory can be confirmed by showing that e 2 is a constant even if the coefficients a and 
b change under different force conditions j41| . The signature of the roughness of the underlying 
energy landscape is uniquely reflected in the non-Arrhenius temperature dependence of the 


unbinding rates. Although it is most straightforward to extract e using Eq.13. no roughness 
measurement, to the best of our knowledge, has been performed using force clamp experiments. 

To extract the roughness scale, e, using dynamic force spectroscopy (DFS) in which the force 
increases gradually in time, an alternative but similar strategy as in force clamp experiments 
can be adopted. A series of dynamic force spectroscopy experiments should be performed as 
a function of T and rf so that reliable unfolding force distributions are obtained. Since a 
straightforward application of EqJT2] is difficult due to the force-dependence of the variables in 
Eq, 12, one should simplify the expression by assuming that the parameters A x*(f), AEq (/), 


and !/£>(/), depend only weakly on /. If this is the case then the second term of Eq, 12 can be 

Arr* 


neglected and Eqjl2] becomes 

k B T 


f* _ k b ± , . k B T 

f K A5‘° g ^ + A? l0E 


+ 


(14) 


' Aad u D e~ AF o/ k BTk B T ' A x$k B T' 

One way of obtaining the roughness scale from experimental data is as follows 021 From the 
/* vs log rf curves at two different temperatures, T\ and T 2 , one can obtain rf{T x ) and r/(T 2 ) 
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for which the /* values are identical. By equating the right-hand side of the expression in Eqfl2l 
at Ti and T 2 the scale e can be estimated mm as 


e 


x 


Axt(Ti)fc B Ti x A x\T 2 )k B T 2 
Axt(Ti)fc B ri - Axt(T2)k B T2 


AF* 




1 1 l tgTl . r/CnjAstpi) 

Axt(T 2 )t i " Axt(T t ) E i^ D (Ti)k B Ti 


k B T 2 
Ax-t (T 2 ) 


log 


r^A^T,) 

VD(T 2 )k B T 2 



Nevo et. al. [42] used DFS to measure e for a biomolecular protein complex consisting 
of nuclear import receptor importin-,/3 (imp-/3) and the Ras-likc GTPase Ran that is loaded 
with non-hydrolyzable GTP analogue (Fig|2]-C). The Ran-imp-/3 complex was immobilized on a 
surface and the unbinding forces were measured using AFM at three values of 77 that varied by 
nearly three orders of magnitude. At high 77 the values of /* increases as T increases. At lower 
loading rates (77 < 2 x 10 3 pN/s), however, /* decreases as T increases (see Fig|2]-B). The data 
over distinct temperatures were used to extract, for the first time, an estimate of e. The values 
of f* at three temperatures (7, 20, 32°G) and EqJT5| were used to obtain e ~ 5 — Qk B T. Nevo 
et. al. explicitly showed that the value of e was nearly the same from the nine pairs of data 
extracted from the f* vs log 77 curves. Interestingly, the estimated value of e is about 0.2 AFq 
where AFg is the major barrier for unbinding of the complex. This shows that for this complex 
the free energy in terms of a one-dimensional coordinate resembles the profile shown in FigjTJ 
It is worth remarking that the location of the transition state decreases from 0.44 nm at 7 °C 
to 0.21 nm at 32 °C. The extracted TS movement using the roughness model is consistent with 
Hammond behavior (see below). 


IV. EXTRACTING TS LOCATION (As*) AND UNFOLDING RATE (k) FROM DY¬ 
NAMIC FORCE SPECTROSCOPY 

The theory of DFS, f* ~ log 77 + log is used to identify the forces that desta¬ 
bilize the bound state of the complex or the folded state of a specific biomolecule. A linear 
regression provides the characteristic extension Ax1 at which the molecule or complex ruptures 
(more precisely Aad is the thermally averaged distance between the bound and the transition 
state along the direction of the applied force). It is tempting to obtain the zero force unfolding 
rate k from the intercept with the abscissa. Substantial errors can, however, arise in the extrap- 
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olated values of Ax* and k to the zero force if /* vs logr/ is not linear, as is often the case when 
rf is varied over several orders of magnitude. Nonlinearity of the [/*, logr/] curve arises for two 
reasons. One is due to the complicated molecular response to the external load that results in 
dramatic variations in Aar. Like other soft matter, the extent of the response (or the elasticity) 
depends on rj |4J S H51 ITT] , The other is due to multiple energy barriers that are encountered in 
the unfolding or unbinding process [31] . 

If the TS ensemble is broadly distributed along the reaction coordinate then the molecule can 
adopt diverse structures along the energy barrier depending on the magnitude of the external 
load. Therefore, mechanical force should grasp the signature of the spectrum of the TS confor¬ 
mations for such a molecule. Mechanical unzipping dynamics of RNA hairpins whose stability 
is determined in terms of the number of intact base pairs is a good example. The conformation 
of RNA hairpins at the barrier top can gradually vary from an almost fully intact structure 
at small forces to an extended structure at large forces. Under these conditions the width of 
the TSE is large. The signature of diverse TS conformations manifests itself as a substantial 
curvature over the broad variations of forces or loading rates. Meanwhile, if the unfolding is 
a highly cooperative all-or-none process characterized by a narrow distribution of the TS, the 
nature of the TS may not change significantly. 

The linear theory of DFS is not reliable if the TSE is plastic because it involves drastic 
approximations of the Eqjl2] From this perspective it is more prudent to fit the the experimental 
unbinding force distributions directly using analytical expressions derived from suitable models. 
If such a procedure can be reliably implemented then the extracted parameters are likely to be 
more accurate. Solving such an inverse problem does require assuming a reduced dimensional 
representation of the underlying energy landscape which cannot be a priori justified. 

V. MECHANICAL RESPONSE OF HARD (BRITTLE) VERSUS SOFT (PLASTIC) 
BIOMOLECULES 

With few exceptions [66], lifetimes of a complex decrease upon application of force. The 
compliance of the molecule is determined by the location of the TS, and hence it is important 
to understand the characteristics of the molecule that determine the TS. As we pointed out, 
many relevant paramters have strong dependence on /, ry , or T. Thus, it is difficult to extract 
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the energy landscape parameters without a suitable model. In this section we illustrate two 
extreme cases of mechanical response da sa im tm of a biomolccule using one-dimensional 
energy profiles. In one example the location of the TS does not move with force whereas in the 
other there is a dramatic movement of the TS. In the presence of force /, a given free energy 
profile F 0 (x ) changes to F(x) = F 0 (x) — fx. The location of the TS at non-zero values of / 
depends on the shape of barrier in the vicinity of the TS. Near the barrier (x « x ts ) we can 
approximate F 0 (x) as 

F 0 (x) « F 0 (x ts ) - (x ts )(x - x ts ) 2 H-• (16) 


In the presence of force the TS location becomes Xt s (f ) = Xt s — y If the transition barrier 
in F 0 (x) is sharp (x fs Fg(xt s ) /) then we expect very little force-induced movement in the 
TS. We refer to molecules that satisfy this criterion as hard or brittle. In the opposite limit 
the molecule is expected to be soft or plastic so that there can be dramatic movements in the 
TS. We illustrate these two cases by numerically computing r^-dependent P(f) using Eq{9]and 
Eqfl2l for two model free energy profiles. 

Hard response: A nearly stationary TS position (independent of /) is realized if the energy 


barrier is sharp (Eq. 16). We model F 0 (x) using 


F 0 (x) = — Vq\(x + l) 2 — £ 2 | with x > 0 


(17) 


where Vo = 28 pN/nm and £ = 4 nm. The energy barrier forms at x = 1 nm and this position 
does not change much even in the presence of force as illustrated in Figj3]4 (top panel). In 
dynamic force spectroscopy the free energy profiles drawn at constant force may be viewed as 
snapshots at different times. The shape of the unbinding force distribution depends on rj. We 
calculated P(f) numerically using Eq|9]and Eqjl2](see the middle panel of Figj3]A.) . Interestingly, 
a plot of the the most probable force /* obtained from P(f) does not exhibit any curvature 
when rj is varied over six orders of magnitude (the bottom panel of Fig|3]A). Over the range 
of rf the [f*,logrf] plot is almost linear. The slight deviation from linearity is due to the 
force-dependent curvature near the bound state (a;&(/)). From the slope we find that Aad ~ 1 
nm which is expected from EqfT7| In addition, we obtained from the intercept in Figj3}-C that 
k = 1.58 s -1 . The value of k[= k(f = 0)] directly computed using Eq|9] is k = 1.49 s -1 . The 
two values agree quite well. Thus, for brittle response the Bell model is expected to be accurate. 
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Soft response: If the position of the TS sensitively moves with force the biomolecule or the 
complex is soft or plastic. To illustrate the behavior of soft molecules we model the free energy 
potential in the absence of force using 

F 0 (x) = —Vo exp (—£x) with x > 0 (18) 

where Vo = 82.8 pN-nm and f = 4 (nm) -1 . The numerically computed P(f) and [/*, log 77 ] plots 
are shown in the middle and bottom panels of Fig|3]B, respectively. The slope of the [/*, logry] 
plot is no longer constant but increases continuously as rj increases. The extrapolated value 
of k to zero / varies greatly depending on the range of ry used. Even in the experimentally 
accessible range of rj there is curvature in the [/*,logrj] plot. Thus, unlike the parameters 
(Ax*, k) in the example of a brittle potential, all the extracted parameters from the force profile 
are strongly dependent on the loading rate. As a result, in soft molecules the extrapolation to 
zero force (or minimum loading rate) is not as meaningful as in hard molecules. Note how the 
extracted Ax* (see the inset in the bottom panel of Fig|3]B) changes as a function of rf. For 
soft (plastic) molecules, the extracted parameters using the tangent at a certain r°j are not the 
characteristics of the free energy profile in the absence of the load, but reflect the features for 
the modified free energy profile tilted by (f*)° at rf. 

In practice, biomolecular systems lie between the two extreme cases (brittle and plastic). In 
many cases the [/*, log rf appears to be linear over a narrow range of ry. The linearity in narrow 
range of log rf, however, does not guarantee the linearity under broad variations of loading rates. 
In order to obtain energy landscape parameters it is important to perform experiments at rf 
as low as possible. The brittle nature of proteins (lack of change in Ax*) inferred from AFM 
experiments may be the result of a relatively large rf (« 1, 000 pN/s). On the other hand, only 
by varying rj over a wide range can the molecular elasticity of proteins and RNA be completely 
described. Indeed, we showed that even in simple RNA hairpins the transition from plastic to 
brittle behavior can be achieved by varying ry [TV]. The load-dependent response may even have 
functional significance. 
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VI. HAMMOND/ANTI-HAMMOND BEHAVIOR UNDER FORCE AND TEMPER¬ 
ATURE VARIATIONS 


The qualitative nature of the TS movement with increasing perturbations can often be an¬ 
ticipated using the Hammond postulate which has been successful in not only analyzing a large 
class of chemical reactions but also in rationalizing the observed behavior in protein and RNA 
foldings. The Hammond postulate states that the nature of TS resembles the least stable species 
along the reaction pathway. In the context of forced-unfolding it implies that the TS location 
should move closer to the native state as / increases. In other words Ax* should decrease as / 
is increased. Originally the Hammond postulate was introduced to explain chemical reactions 
involving small organic molecules EH E3J. Its validity in biomolecular folding is not obvious 
because there are multiple folding or unfolding pathways. As a result there is a large entropic 
component to the folding reaction. Surprisingly, many folding processes are apparently in accord 
with the Hammond postulate [ 781180 ] . If the extension is an appropriate reaction coordinate for 
forced unfolding then deviations from Hammond postulate should be an exception rather than 
the rule. Indeed, anti-hammond behavior (movement of the TS closer to more stable unfolded 
state as T increases) was suggested by Schlierf and Rief [43] based on a model used to analyze 
the AFM data. The simple free energy profiles used in the previous section (Eq jTr] and Eq, 18) 
can be used to verify the Hammond postulate when the external perturbation is either force or 
temperature. First, for the case of hard response the TS is barely affected by force, thus the 
Hammond or anti-Hammond behavior is not a relevant issue when unbinding is induced by /. 
On the other hand, for the case of soft molecules Ax* always decreases with a larger force. The 
positive curvature in [/*,log7y] plot is the signature of the classical Hammond-behavior with 
respect to /. 

As long as a one dimensional free energy profile suffices in describing forced-unfolding of 
proteins and RNA the TS location must satisfy the Hammond postulate. In general, for a fixed 
force or ry, Ax* can vary with T. The changes in Ax* with temperature can be modeled using 
T-dependent parameters in the potential. To evaluate the consequence of T-variations we set 


f = £ 0 + a(T - 300A') (19) 

for both free energies in Eq|l7]and Eq]l8] Depending on the value of a the position of the TS can 
move towards or away from the native state. We set a = ±0.1 for both the hard and soft cases. 
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The numerically computed [/*,logrj] plots are shown in Figj3} One interesting point is found 
in the soft molecule that exhibits Hammond behavior. For wide range of 77 , Aad decreases as 
T increases. However, the most probable unbinding force f* at low temperatures can be larger 
or smaller than f* at high temperatures depending on the loading rate (see upper-right corner 
of Figj3]). A very similar behavior has been observed in the forced-unbinding of Ran-imp-/3 
complex [42] (see also Figj2]-B). Although the model free energy profiles can produce a wide 
range of behavior depending on T, /, and 77 the challenge is to provide a structural basis for 
the measurements on biomolecules. 

VII. MOLECULAR TENSEGRITY AND THE TRANSITION STATE. 

For Rts (see FigjijA for the definitions of Rts, Ru, Rf on a one-dimensional free energy 
profile F(R)), associated with the barrier top of F(R) at f — f m where f m is the transition mid¬ 
force, to be considered the “true” transition state, it is necessary to ensure that it is consistent 
with other conventional definitions of the transition state ensemble. A plausible definition of 
the TS is that the forward (to the unfolded state, P un foid ) and backward (to the folded state, 
Pfoid ) fluxes starting from the transition state on the reaction coordinate should be identical 
[ 8 TJ . For an RNA hairpin it means that if an ensemble of structures were created starting 
at Rts then the dynamics in the full multidimensional space would result in these structures 
reaching the folded and unfolded states with equal probability. The number of events reaching 
Rf and Ru starting from Rts can be directly counted if folding trajectories with high temporal 
resolution exhibiting multiple folding and unfolding transitions at / = f m can be generated 
(FigJdjA.). Our coarse-grained simulations of the P5GA RNA hairpin, which were the first to 
assess the goodness of Rts as a descriptor of the TS, showed that starting from Rts the hairpin 
crosses the TS region multiple times before reaching R = R F or R F , suggesting that the TS 
region is broad and heterogeneous. The transition dynamics of biopolymers occurs on a bumpy 
folding landscape with fine structure even in the TS region, which implies there is an internal 
coordinate determining the fate of trajectory projected onto the A-coordinate. In accord with 
this inference we showed that for the P5GA hairpin the TS structural ensemble is heterogeneous 
(Fig. [4^3). More pertinently, the forward and backward fluxes starting from the structure in TS 
ensemble (see the dynamics of trajectories starting from Rts in Fig|4)3) do not satisfy the equal 
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flux condition, Pfoid = Punfoid = 0.5. Thus, from a strict perspective R for a simple hairpin 
may not be a good reaction coordinate even under tension, implying that R is unlikely to be an 
appropriate reaction coordinate at / ~ 0. 

Based on simulations Morrison et al. [82] proposed a fairly general theoretical criterion to 
determine if R could be a suitable reaction coordinate. The theory uses the concept of tensegrity 
(tensional integrity), which was introduced by Fuller and developed in the context of biology to 
describe the stability of networks. The notion of tensegrity has been used to account for cellular 
structures [83] and more recently for the stability of globular proteins [84), the latter of which 
made an interesting estimate that the magnitude of iuter-residue precompression and pretension, 
associated with structural integrity, can be as large as a few 100 pN. Using F(R) the experimen¬ 
tally measurable molecular tensegrity parameter is defined s = f c /f m = AiA(/ m )// m AiU(/ m ), 
where A F* = F(Rts) — F(Rp) and A R^ = Rts — Rf • The molecular tensegrity parameter s 
represents a balance between the compression force (/ m ) and the tensile force (/ c ), a building 
principle in tensegrity systems [85]. For hairpins such stabilizing interactions involve favor¬ 
able base pair formation. In terms of s and the parameters characterizing the one dimensional 
landscape (/ m and k u in Fig. [4jA) an analytic expression for P un f 0 id has been obtained [S2J. 


Punfold(s) 


1 1 
2 1+s 

1 r $ i 32/,u$—i) ■ 

2 l+$ “ r nk u (<f>+l) 3 * 


S > 1 
S < 1 


( 20 ) 


For R to be a good reaction coordinate, it is required that P un foid ~ The theory has 
been applied to hairpins and multi-state proteins. Using experimentally determinable values 
for s, one can assess whether R is a good reaction coordinate by calculating P un foid using 
theory. Applications of the theory to DNA hairpins [82] showed (Fig. |4p) that the precise 
sequence determines whether R can be reliably used as an appropriate reaction coordinate, thus 
establishing the usefulness of the molecular tensegrity parameter. Application of the theory 
based on molecular tensegrity showed that R is not a good reaction coordinate for riboswitches 

m- 
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VIII. MULTIDIMENSIONALITY OF ENERGY LANDSCAPE COUPLED TO 
’’MEMORY” AFFECTS THE FORCE DYNAMICS 

The natural one dimensional reaction coordinate in mechanical unfolding experiments is the 
extension x of the molecule. In studies of chemistry and biophysics, it is quite standard to 
project the dynamics of a molecule with many degrees of freedom onto a one-dimensional (ID) 
reaction coordinate. Especially, in single molecule force experiments, the molecular extension 
“x”, parallel to the direction of an external force, is routinely employed to describe the force- 
induced dynamics of biomolecules. However, the assumption of ID projection is valid only if 
time scales of dynamics is clearly separated between a variable associated with the reaction 
coordinate and other degrees of freedom [HS]. For a complex biomolecule that has multiple 
time scales intertwined in its dynamics, the signature of interference between two variables 
with comparable time-scales should be displayed in experimental measurements. As a simple 
extension of the theories discussed above, we consider a situation of an unbinding process where 
the fluctuation of a molecular surface between open and closed states can gate the unbinding 
kinetics of a ligand from its binding pocket (Fig.2). Depending on the ratio kj A where A is 
the rate of gating and k is the time scale for unbinding in the absence of gating, the ligand is 
expected to undergo disparate unbinding kinetics. If k/X S> 1 or k/X <C 1, the environment 
appears static to the ligand |87, 88] - Thus, the ligand unbinding occurs via parallel paths over 
multiple barriers (k/X 1) or via single path over a rapidly averaging barrier (k/X <C 1). 
Whereas, if k/X ~ 0( 1), the open-H-closed gating produces a fluctuating environment along the 
dynamic pathway of the ligand and affects the unbinding process in a non-trivial fashion, which 
is often termed dynamical disorder. Suppose that the reaction rate constant associated with 
barrier crossing is controlled by the openness of the receptor molecule, which can be modeled 
as kr 2 where r is the radius of the bottleneck. Provided that the distribution of bottleneck size 
is Gaussian, i.e., <p(r) ~ e ~ r2 / 26 with 9 = (r 2 ), and r obeys a stochastic differential equation 
d t r = — Ar+£ r (t) with (£ r (i)£ r (i / )) = 2A 95(t—t'), two disparate unbinding kinetics are quantified. 
If A —> oo much greater than barrier crossing rate then the reactivity is defined by a pre-averaged 
reaction constant with respect to r, i.e., kr 2 —y kd. In this case, the survival probability of 
the ligand decays exponentially (S(t) = e ~ k6t ) with k = koe^ xt . In contrast, if A —» 0 
(A -C k ) so that the memory of the initial r value is retained while barrier crossing takes place, 
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S(t) = J 0 °° e~ kr2t (p(r)dr ~ (1 + 2k9t)~ 1 ^ 2 . Two scenarios of bottleneck gating frequency lead 
to totally different consequences. The gating frequency A, intrinsic to a molecule of interest, is 
difficult to adjust independently although it can in principle be varied to a certain extent by 
changing viscosity. However, in single molecule force spectroscopy it is feasible to vary k by 
controlling the external tension, so that the ratio k/\ can be scanned from fc/A^ltofc/A<^l 
via k /A ~ 0(1). 

The physical picture described above can be mathematically formulated as follows. In the 
context of ligand binding to myoglobin Zwanzig first proposed such a model by assuming that 
reaction (binding) takes place along the ^-coordinate and at the barrier top (x = x ts ) the 
reactivity is determined by the cross section of bottleneck described by the r-coordinate [fSTj . 
We adopt a similar picture to describe the modifications when such a process is driven by force. 
First, consider the Zwanzig case i.e, with loading rate rf — 0. The equations of motion for x 
and r are, respectively, 


m 


d 2 x 

dt 2 

dr 

dt 


dx 

dt 


dU(x ) 


dx 

—7 r + F r(t). 


F x (t ) 


( 21 ) 


The Liouvillc theorem (^ = 0) describes the time evolution of probability density, p(x,r,t ), as 


| = f + ^fe) + ^(§P]=°. 

dt dt dx \dt ) or \dt 


( 22 ) 


By inserting of Eq.21 to Eq.22 and neglecting the inertial term, (m^f), and averaging over the 
white-noise spectrum, and the fluctuation-dissipation theorem ((F x (t)F x (t')) = 2 (ksTS(t — t'), 
(F r (t)F r (t')) = 2A 95(t — t') where (r 2 ) = 9) leads to a Smoluchowski equation for p(x,r,t ) in 
the presence of a reaction sink. 


^ = C x p + C r p - k r r 2 5(x - x ts )p 


(23) 


where C x = + and C r = A 9-^ + |). Integrating both sides of Eq 

using J dxp(x,r,t ) = C(r,t ) leads to 
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d(J 

~dt 


C r C(r , t) - k r r 2 p(x ts ,r , t). 


(24) 


By writing p(x ts ,r,t ) = 

e -U(x ta )/k B T j J ^ xe ~U(x)/k B T 


( f> x (xts)C(r,t ) where 4>(xt s ) should be constant as 4>(xt s ) = 

becomes 


U"(x b ) r -(U(x ta )-U(x h ))/k B T 


27rfc^T 


, Eq 
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dC 

~dt 


C r C(r,t) — kr 2 C(r,t). 


(25) 
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where k = k r \J^-^^ e -( u ( x ti>)-u(x b ))/k B T _ j n a jj likelihood k r reflects the dynamics near the 
barrier, so we can write k = k e~ AU / k B T where n describes the geometrical information of the 
cross section of bottleneck. Now we rederive the equation in Zwanzig’s seminal paper where the 
survival probability (E (t) = J 0 °° drC(r , £)) is given under a reflecting boundary condition at r = 0 
and Gaussian initial condition C(r,t = 0) ~ e~ r2 1 26 . By setting C(r,t ) = exp (v(t) — n(t)r 2 ), 


Eq.25 can be solved exactly, leading to 


u'{t) = —2X6 n{t) + A 

n'(t) = —AX6fi 2 {t)-\-2Xfi{t)-\-k. (26) 

The solution for /i{t) is obtained by solving y Jl'lfo ' a 2 _iQ 9 2 a 2 = t, and this leads to 

= i r js + i)-(s-i)e 

MO) 2 \ (S + 1) + (S-1)E 

, w = M ( g - 1)+ ,° g ( (5+1) tf~ 1)g )~ 1/2 < 27 > 

with /i(0) = 1/20. The survival probability, which was originally obtained by Zwanzig, is 

E(t) ^ xp (^ (5 - 1)f )( (^^-/-^ )- 1/2 (28) 

where S = (l + ^) V2 and E = e~ 2XSt . 

We wish to examine the consequences of coupling between local and global reaction coordi¬ 
nates under tension. In order to accomplish our goal we solve the Smoluchoski equation in the 
presence of constant loading rate. In this case, k in Eqj25]should be replaced with f te t ( r f / ^ xt /k B T)^ 


Eq.25, however, becomes hard to solve if the sink term depends on t. Nevertheless, analytical 
solutions can be obtained for special cases of A. If A —> oo, dE(t)/dt = —kde t ^ r f Axt ^ kBT ^ E(t), 
and hence, 


E(/) = exp 


-mJ o = exp _ i)^ (29) 


Using the rupture force distribution P(f) = — and dP jp \f=f* = 0, one can obtain the most 
probable force 

kr, T TrA'rt 

(30) 


r EEL log r ’ Axt 


Ax t & ( kd)k B T' 
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If A is small (A —* 0) then C(r,t ) = exp [— dtkr 2 exp (t x rfA/k B T)\ = 


exp 


kr 2 


ryAtct 


|«T( e tr / Axt/fcBr _ ^ with the initial distribution of e"” 2 / 26 . Thus, 


S(/) = / dr exp 



r 2 kk B T , c fAx*/k R T _ 1 ' ) 

rfAxi 

( kkjff , fAx t /taT __.j, 1t 1/2 


- exp [-r 2 /29] 


(31) 


Note that if ry —* 0 we recover Zwanzig’s result £(i) ~ (1 + 2A;6^) 1 / 2 . Using P(f) = —dT,(f)/df 


p( n = 1 f kkBT ( r f^/k B T _ N , JA 3/2 kfAx*/k B T 

[j) s/2^e\r f Axt [e ) + 2 e) r f e 


dP(f)/df\f =f * = 0 gives 


f* = 


k bT ( / rfAx ^ A 


Ax$ 


log 


2 kdksT 


\kdk B Tj V r f Ax* 


(32) 


(33) 


m 


which rf > (1 + 26) k9 ^t" since /* > 0. This shows that /* vs rj has a different form when 


A —» 0 from the one when A —> oo. The deviation of Eq.33 from the conventional relation is 
pronounced when [77 — (1 + 26) —> 0 + . 

At present, experimental data have been interpreted using mostly one dimensional free en¬ 
ergy profiles. The meaning and the validity of the extracted free energy profiles has not been 
established. At the least, this would require computing the force-dependent first passage times 
using the ’’experimental” free energy profile assuming that the extension is the only slowly re¬ 
laxing variable. If the computed force-dependent rates (inversely proportional to first passage 
times) agree with the measured rates then the use of extension of the reaction coordinate would 
be justified. In the absence of good agreement with experiments other models, such as the one 
we have proposed here, must be considered. In the context of force-quench refolding we have 
shown (see below) that extension alone is not an adequate reaction coordinate [89j. For refolding 
upon force-quench of RNA hairpins, the coupling between extension and local dihedral angles, 
which reports on the conformation of the RNA, needs to be taken into account to quantitatively 
describe the refolding rates. 


IX. CONCLUSIONS 


With the advent of single molecule experiments that can manipulate biomolecules using 
mechanical force it has become possible to characterize energy landscapes quantitatively. Me- 
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chanical folding and unfolding trajectories of proteins and RNA show that there is great diversity 
in the explored routes 12311351 190 ). In certain well defined systems with simple native states, 
such as RNA and DNA hairpins, it has been shown using constant force unfolding that the 
hairpins undergo sharp bistable transitions from folded to unfolded states [23, 23 08J- From 
the dynamics of the extension as a function of time measured over a long period the underlying 
force dependent profiles have been inferred. The force-dependent folding and unfolding rates 
and the unfolding trajectories can be used to construct the one-dimensional energy landscape. 
In a remarkable paper [25], Block and coworkers have shown that the location of the TS can be 
moved, at will, by varying the hairpin sequence. The TS was obtained using the Bell model by 
assuming that the Ax* is independent of /. While this seems reasonable given the sharpness of 
the inferred free energy profiles near the barrier top it will be necessary to show the Ax* does 
not depend on force. 

The fundamental assumption in inverting the force-clamp data is that the molecular extension 
is a suitable reaction coordinate. This may indeed be the case for force-spectroscopy in which 
the response of the molecule only depends on force that is coupled to the molecular extension, 
which may well represent the slow degrees of freedom. The approximation is more reasonable 
for forced-unbinding. It is less clear if it can be assumed that extension x is the appropriate 
reaction coordinate when refolding is initiated by quenching the force to low enough values 
such that the folded state is preferentially populated. In this case the dynamic reduction in 
x can be coupled to collective internal degrees of freedom. In a recent paper m we showed, 
in the context of force-quench refolding of an RNA hairpin, that the reduction in x is largely 
determined by local conformational changes in the dihedral angle degrees in the loop region. 
Zipping by nucleation of the hairpin with concomitant reduction in x does not occur until the 
transitions from trans to gauche state in a few of the loop dihedral angles take place. In this 
case, one has to consider at least a two dimensional free energy landscape. Figj 6 ] clearly shows 
such a coupling between end-to-end distance (R) and the dihedral angle degrees of freedom. The 
’’correctness” of the six dihedral angles representing the conformation of the RNA hairpin loop 
region (0j, i — 19,... 24) is quantified using (1 — cos (0* — 0°)), where 0° is the angle value in the 
native state and (...) is the average over the six dihedral angles. (1 — cos (0 — 0°)) = 0 signifies 
the correct dihedral conformation for the hairpin loop region. Once the ’’correct” conformation 
is attained in the loop region, the rest of the zipping process can easily proceed as we have 


23 


shown in m- Before the correct loop conformation is attained, RNA spends substantial time 
in searching the conformational space related to the dihedral angle degree of freedom. The 
energy landscape ruggedness is manifested as in Fig|6]when conformational space is represented 
using multidimensional order parameters. The proposed coupling between the local dihedral 
angle degrees of freedom and extension (global parameter) is fairly general. A similar structural 
slowing down, due to the cooperative link between local and global coordinates, should be 
observed in force-quench refolding of proteins as well. 

One of the most exciting uses of singe molecule experiments is their ability to extract precise 
values of the energy landscape roughness e by using temperature as a variable in addition to 
/. In this case a straightforward measurement of the unbinding rates as a function of / or 
rf can be used to obtain e without having to make any assumptions about the underlying 
mechanisms of unbinding. Of course, this involves performing a number of experiments. In so 
doing one can also be rewarded with a diagram of states in terms of / and T m- The theoretical 
calculations and arguments given here also show that the power of single molecule experiments 
can be fully realized only by using data from the experiments in conjunction with carefully 
designed theoretical and computational models [5TJ. The latter can provide the structures that 
are sampled in the process of forced-unfolding and force-quench refolding as was illustrated for 
ribozymes and GFP [90]. It is likely that the promise of measuring the energy landscapes of 
biomolecules, one molecule at a time, will be fully realized using a combination of single molecule 
measurements, theory, and simulations. Recent studies have already given us a glimpse of that 
promise with more to come shortly. 
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Figure Caption 


Figure [l] : Caricature of the rough energy landscape of proteins and RNA that fold in an 
apparent ”two-state” manner using extension x, the coordinate that is conjugate to force /. 
Under force /, the zero-force free energy profile (F 0 (x)) is tilted by / x x and gives rise to the 
free energy profile, F(x ). In order to clarify the derivation of Eq(7]we have explicitly indicated 
the average location of the relevant parameters. 

Figure [2] : Dynamic force spectroscopy measurements of single imp-/3- RanGppNIIp pairs 
at different temperatures. A. Distributions of measured unbinding forces using AFM for the 
lower-strength conformation of the complex at different loading rates at 7 and 32 °C. Roughness 
acts to increase the separation between the distributions recorded at different temperatures. 
The histograms are fit using Gaussian distributions. The width of the bins represents the 
thermal noise of the cantilever. B. Force spectra used in the analysis. The most probable 
unbinding forces /* are plotted as a function of log(r/). The maximal error is ±10% because of 
uncertainities in determining the spring constant of the cantilevers. Statistical significance of 
the differences between the slopes of the spectra was confirmed using covariance test. (Images 
courtesy of Reinat Nevo and Ziv Reich [45]). C. Ran-importin/3 complex crystal structures 
(PDB id: 1IBR [92j) in surface (left) and ribbon (right) representations. In AFM experiments, 
Ran (red) protein complexed to irnportin/3 (yellow) is pulled until the dissociation of the 
complex takes place. 

Figure [3] : Dynamic force spectroscopy analysis using two different free energy models. 
From the top to bottom, (top) free energy profiles at different forces, (middle) rupture force 
distributions, and (bottom) relationship of most probable force as a function of log-loading 
rates are shown for (A) hard and (B) soft potentials. For hard potential, free energy profile 
in the absence of force is F 0 (x) = —Vq\(x ± l) 2 — £ 2 j with Vb = 20 pN/nm, £ = 4 nm, and 
x > 0. The lack of change in x ts as / changes shows a brittle response under tension. For soft 
potential, free energy profile in the absence of force is F 0 (x) = — Vbexp(— £x) with V 0 = 82.8 
pN-nrn, £ = 4 (nm) -1 . The lack of change in xt s as / changes shows a brittle response under 
tension. For emphasis on the soft response of the potential, the position of TS at each force 
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value is indicated with arrows. C. By tuning the value of £ (see Eq. (23)) as a function of 
temperature, Hammond and anti-Hammond behaviors emerge in the context of force spectra 
in the free energy profiles that show hard and soft responses. The condition for Hammond or 
anti-Hammond behavior depends on a (Eq. (23)). 

Figure [5} Fluctuating bottleneck model under an external force. Gating dynamics of 
receptor protein conformation with frequency A can interfere with the dynamics of ligand 
unbinding if the unbinding rate (k) is comparable or slower than A. 

Figure [4j A. Time trace of end-to-end distance (R) of RNA hairpin at transition midforce 
f m = 14.7 pN (left). The corresponding free energy profile in terms of R, F(R) (right). The 
positions of native, unfolded, and transition states are marked with arrows. In addition, barrier 
height (A F*) and the curvature of the unfolded state (ku) are also shown on the F(R). B. Time 
trajectories of simulations starting from the configurations of the transition state ensemble 
(shown on the left). Trajectories reaching the folded and unfolded state at 2.0 and 7.5 nm are 
colored in blue and red, respectively. In the blue trajectories, a number of recrossing events 
can be observed. C. Tensegrity parameters calculated for four DNA hairpins with different 
sequences in Ref. [93] are related to P U nfoid■ The DNA hairpin with sequence B is predicted 
to have Punfoid most proximal to 0.5, which suggests that the free energy profile calculated in 
terms of end-to-end distance coordinate most accurately describes the dynamics of this DNA 
hairpin. 

Figure [6] : A. A sample refolding trajectory of a RNA hairpin starting from the stretched 
state. The hairpin was, at first, mechanically unfolded to a fully stretched state and the force 
was subsequently quenched to zero at t ~ 20 {is. The time-dependence of the end-to-end 
distance shows that force-quench refolding occurs in steps. B. The deviation of the dihedral 
angles from their values in the native state as a function of time shows large departures from 
native values of the dihedral angles in loop region (indicated by the red strip). Note that this 
strip disappears around t ~ 300 /is, which coincides with the formation of bonds shown in C. fs 
is the fraction of bonds in pink that indicates that the bond is fully formed. D. The histograms 
collected from the projections of twelve stretching and force-quench refolding trajectories on the 
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two dimensional plane characterized by the end-to-end distance (R) and the average correctness 
of dihedral angles ((1 — cos (0 — 4> 0 ))) around the loop region (i — 19:*— 24). The scale on the 
right gives the density of points in the two dimensional projection. This panel shows that the 
local dihedral angles are coupled to the end-to-end distance R , and hence extension alone is not 
a good reaction coordinate especially in force-quench refolding. The molecular extension x is 
related to R by x = R — Rn where Rn is the distance in the folded state. 
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